Shear-free gravitational collapse is strongly censored 
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We consider spherically symmetric spacetimes with matter whose timelike flow is assumed to be 
shear- free. A number of results on the formation and visibility of spacetime singularities is proven, 
with the main one being that shear-free collapse cannot admit locally naked singularities (which 
implies absence of globally naked singularities). We conjecture that shear is a necessary condition 
for the occurrence of locally naked singularities in generic gravitational collapse. 
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I. INTRODUCTION 



The final state of gravitational collapse is one of the 
outstanding issues in classical general relativity. One par- 
ticular aspect of this problem is the role of shear, or 
absence thereof, in the formation and visibility of sin- 
gularities. Here, we investigate the formation and local 
visibility of singularities in spherical shear-free collapse. 
Our motivation is fourfold: (i) Gravitational thermody- 
namics considerations suggest that shear affects the end- 
state of collapse [l]; (ii) Spherical shear-free spacetimes 
constitute a vast class, which include cosmological (i.e., 
non-asymptotically flat) solutions Q, isolated soliton-like 
objects [3, and isolated collapsing objects [I]; (iii) Rel- 
ativistic kinetic theory tells us that, when the distribu- 
tion function is isotropic, the solution to the Einstein- 
Boltzmann equations must be shear- free Q; (iv) Shear- 
free (perfect fluid) solutions constitute a useful model for 
astrophysical collapse of globular clusters [||. 

Within spherical symmetry, studies of certain shear- 
free collapse models 0, 0, H, EJ Ell > as well as Tolman- 
Bondi (nonvanishing shear) collapse ^U, suggest that, 
provided certain physical reasonability conditions are sat- 
isfied, when the matter flow is shear-free resulting space- 
time singularities are not naked. In this paper, we 
consider spherical spacetimes with a general anisotropic 
fluid, whose timelike flow is assumed to be shear-free; the 
present analysis includes and generalizes all the spherical 
shear-free collapse models in the literature thus far. 

The main purpose of this paper is threefold: (i) De- 
rive explicit analytical conditions for geodesic complet- 
ness of the spacetime; (ii) Show that shear-free collapse 
always leads to a spacetime singularity; (iii) Show that 
such sing ularity is locally covered — and thus globally cov- 
ered — along all causal (radial and nonradial) direc- 
tions. The latter result is the main one in the paper, and 
may be viewed as a strong cosmic censorship theorem 
for spherical shear- free spacetimes. We also discuss the 
interplay between shear and Weyl curvature, and put for- 
ward the conjecture that shear is a necessary condition 
for singularities to be (at least) locally visible in generic 
gravitational collapse. By "generic" , we mean a space- 
time of arbitrary geometry, which admits a well-posed 
regular Cauchy data problem, and whose matter content 



obeys the weak energy condition. Natural geometrizcd 
units, in which 8nG = c = 1, are used throughout. 



II. SPHERICAL SHEAR-FREE SPACETIMES 

We start with the most general spherically symmetric 
metric: 



ds 1 



2>1< 



(It 2 



z^dr 2 



R 2 {d8 



2 ' sin 2 



! ), (2.1) 



where t G R, r G 9 G [0,tt], ip G [0,2tt), and $, 
and R are functions of t, r alone. For the matter content 
we take the stress-energy tensor: 



T£ = diag(-p,n,r,r) ) 



(2.2) 



where II and T are the principal pressures, parallel and 
orthogonal to the timclikc fluid flow = e~*<5f , respec- 
tively. The field equations may be written as: 



R,tr 

n r 

M(t,r) 
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1 
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2(r-n)^- (/ . + niM/.,. 



(2.3) 
(2.4) 

(2.5) 

(2.6) 

(2.7) 



where the last equation is the Misner-Sharp mass |13j . 
Solutions with the stress-energy tensor l|2.2|l are fully 
characterized by the kinematical properties of the 4- 
velocity u^: vorticity w M „, acceleration a M , expansion O, 
and shear a^: 



Uuv = 0, 



e = 



($ t + 2R.tR- 1 ), 
of = diag(0,-2er,cr,cr), 



(2.8) 
(2.9) 
(2.10) 
(2.11) 



where the vorticity vanishes because of spherical symme- 
try, and 



1 



* = -e-^ORtiT 1 -^). 



(2.12) 
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Shear-free motion is given by R.tR 
grates to 

e* = f(r)R, 



where f(r) is an arbitrary positive-definite real- valued 
function, which we set as f(r) = l/r without loss of 
generality. Equation 1|2.3[1 then integrates to 



2* 



(R t ) 2 

R 2 H 2 (t)' 



(2.14) 



where H{t) is (one third of) the expansion of the flow 
lines, as can be seen from Eq. H2.10|) . For the case of per- 
fect fluids (n = r) , the remaining field equations lead to 
a single second-order nonlinear partial differential equa- 
tion Q: 



F xx + J{x)F 2 = 0, 



where x := r 2 , F(t,r) := yfx/R, and J(x) is an arbi- 
trary function. Equation (|2.15|l fully determines the so- 
lution to spherical shear-free perfect fluid collapse. For 
J.x = 0, the PDE can be solved in terms of Weierstrass 
p- functions 0, and physically it corresponds to the so- 
called synchronous evolution, wherein all shells collapse 
to vanishing proper area in the same amount of proper 
time 0. For J iX ^ 0, Eq. I|2.15|l cannot in general be 
solved in closed form, and very few cases admit analytical 
solutions . In this paper we shall consider the general 
case (IT ^ r), noting, however, that all of the subsequent 
results also hold for the perfect fluid case. 

We now impose the following regularity conditions on 
the metric: 

R(t,0) = i? t (i,0) = 0, (2.16) 
limr— = 1, (2.17) 
3 AlI/ eC 1 (R) and |5 wa | r =o < oo. (2.18) 

We note that, although C 2 differentiability is required 
for non-singular Riemann curvature, we shall not impose 
it (other than on the initial Cauchy data, as detailed 
below) , since one would like to allow for curvature singu- 
larities, which arise dynamically from the gravitational 
collapse of regular Cauchy data; in addition, local flat- 
ness and geodesic motion only require g^ E C 1 (R), since 
only the Christoffcl connection enters the geodesic equa- 
tion. 

In spherical symmetry there are no gravitational de- 
grees of freedom, and the metric is therefore uniquely 
determined from Cauchy data on a spacelike hypersur- 
face T, tg — {x^ : x° = to = const.}. To ensure that 
any putative singularities arise dynamically from regular 
initial data, we impose the following conditions on it: 



$ t, which inte- Condition (|2.19() ensures smooth Riemann curvature on 
S to , and (|2.20(l guarantees the initial slice is free of 
trapped surfaces. 
(2.13) For the collapse of isolated bodies, the metric i|2.1|) 

must be matched to a Schwarzschild exterior at some fi- 
nite coordinate radius r c , via the Darmois-Israel junc- 
tions conditions, which impose continuity of the met- 
ric and extrinsic curvature across the junction surface. 
Since we shall be interested in the R — singularity aris- 
ing from gravitational collapse of any spherical shear-free 
spacetime — not necessarily asymptotically flat — we shall 
not discuss the details of the matching (see, e.g., [l7|l. 
which fixes the function H (t) for a given interior solution 
$, R, p, IT, r}. For collapsing solutions R.t < and 
H(t) < 0, and the particular functional form of H (t) is 
qualitatively irrelevant for the endstate of collapse. 

Before proceeding, we compute the Brown- York 
quasilocal mass |18fl . which provides a useful tool to an- 
(2.15) alyze the role of "energy" in the formation of curvature 
singularities. Consider a spacelike hypersurface S t with 
two-boundary dT, t ; the quasilocal mass E on E t associ- 
ated with the three-volume defined by dH t is the value 
of the Hamiltonian that generates unit time translations 
orthogonal to the boundary dT, t [Isj |. In terms of the in- 
duced two-metric j a b on <9£ t , it is given by the proper 
surface integral 



E:=±- 
8tt 



d 2 x^{K - K ), 



(2.21) 



where K is the trace of the extrinsic curvature of 9E t as 
imbedded in St, and Kq is the corresponding quantity 
for the same surface as embedded in a three-dimensional 
slice of flat spacetime (such embedding is assumed to be 
possible). This last term corresponds to a normalization 
of the "zero of energy" with respect to flat space. 
In our case, the Riemannian metric on E t is 



hijdx l dx J = e 



■24' 



dr 2 + R 2 d9 2 



R 2 sin 2 



(2.22) 



For the two-boundary <9£ t we take the surface S, para- 
metrically defined by E(x^) = r — ro = 0, which has the 
induced Riemannian two-metric: 

lab dw a dw b = R 2 {t,r a )d9 2 + R 2 (t,r ) sin 2 9 dip 2 . (2.23) 

The extrinsic curvature of S as embedded in T, t is 



K ab = n. 



( d 2 



\dw a dw b 



_i_ ( 3 )r i 



dxi dx k 



Jk dw a dw b 



where 



n//i ij ' (3) ViH( 3 )Vj-H 



3*ff 



(2.24) 



(2.25) 



g^(h,r) eC 2 (R), 
2M(t ,r) < R(t ,r). 



(2.19) 
(2.20) 



is the unit normal to S. The only non-vanishing compo- 
nents of K a b are 



K vv = sin 2 9K ee = —e~®RR r , 



(2.26) 



3 



whence 



K = K abl ab 



-2e 



_<!> R,r 

~R 



(2.27) 



For the Euclidean reference surface, we take the metric 

~ 2 (2.28) 



E ij dx i dx j = dr 2 + r 2 d9 2 



r 2 sin 2 



A trivial calculation gives K = — 2/r, and the Brown- 
York mass is then 



£ = 



1 

8^ 



dip / d9r 2 sinO -2e 



r 1 — e r— — 



(2.29) 



where all quantities are evaluated at r = ro on E t . This 
quasilocal mass agrees with the ADM and Bondi-Sachs 
mass in the appropriate limits, and a "positive-mass" 
theorem for it has recently been proven [19j . To make 
contact with the Schwarzschild spacetime (with mass pa- 
rameter m) , set R — r and <f> = — In \/T— 2m/ r in Eq. 
(|2.29(l . whereby Eq. (6.14) in 18] is recovered. This 
shows that E asymptotes m (and thus M) in the weak- 
field limit, m/r <C 1. 



III. CAUSAL GEODESIC COMPLETENESS 

A spacetime is said to be causally g-complete if ev- 
ery non-spacelike geodesic 7(A) can be extended to arbi- 
trary values of its affine parameter |20| . One must note 
that g-completeness does not imply that the spacetime 
is singularity-free, since non-geodesic (but causal) inex- 
tendible curves can exist in g-complete spacetimes |2l|. 
A spacetime is singularity-free if all C 1 curves, geodesic 
or otherwise, are extendible, in which case it is said to 
be b-complete; b-completeness implies g-completeness, 
but the converse is not true. Evidently, then, failure of 
g-completeness is sufficient to show that the spacetime 
contains a singularity. From elementary ODE theory, a 
necessary and sufficient condition for g-completeness is 
the boundedness of all of the first derivatives of the coor- 
dinates with respect to the affine parameter j^]. Since 
we are interested in showing that shear-free collapse al- 
ways leads to a spacetime singularity, it suffices to take 
a particular family of causal geodesies and show that it 
is inextendible, i.e., that at least one of the functions 
dx 11 /dX grows unbounded for finite A. 

We start with a generic geodesic, 7(A), defined by the 
integral curves of the vector field 



dr* 1 



(3.1) 



where A is an affine parameter along 7(A), and the over- 
dot denotes total differentiation with respect to A. Nor- 
malizing = — e, we have 



i 2 = e" 2 * 



2* -2 

e r 



L 2 
R 2 



(3.2) 



where L is the total conserved angular momentum, given 
by L 2 = L 2 e) +L 2 v) , where 



J (6) - °e 



5%^ = 6R 2 , L (v) = Sit^ = ipR 2 sin 



(3.3) 



We remark that, without loss of generality, one can al- 
ways take the polar coordinates to be such that 7(A) 
moves along the 8 — tt/2 hypersurface, whence L = ipR 2 ; 
obviously, this does not change Eq. (|3.2(l . nor the 
geodesic equation: 



-2* 



R 



(3.4) 



2t + t 2 (* f + $ ( ) + 2if* r + — e 

R z 

- e e- 2 *$ if = 0, 
2r + f 2 (* r + $, r ) + 2tr$ t -!J e - 2 * (^f " *,r) 

- ee" 2 ** r = 0. (3.5) 

Equations l|3.2|l . 1)3.411 - 1)3.5(1 form a complete set for the 
three variables {t, r,(p}. 

We now specialize to the case of future-oriented ingoing 
radial null geodesies (IRNG): 

i> 0, r < 0, 6 = <f = e = 0, (3.6) 
2i'+t 2 (1', t + $ it ) + 2ir*,r = 0, (3.7) 
2f + r 2 (# r + $r) + 2tr$ t = 0. (3.8) 

The null condition e = reduces these two equations to 
a single first-order ODE for t: 



1 di 

i dX ~ 

which integrates to 



* + $ +*($,< - 



t = (3(r)e- 



(3.9) 



(3.10) 



where (3(r) : > R + \{0} is an arbitrary integration 

function. We can now put forward the following results: 
Proposition 1. A spherically symmetric spacetime with 
metric \2.1)) is g-complete along geodesies defined by 
\3.6}) - \3.§) . if the following holds: 

3C e M+\{0} : r) > lnC + /(r), Vi e R, Vr e R+ 

Proof. The condition above implies that t remains fi- 
nite and bounded away from zero. Since \r\ = ie*~*, 
and e C 1 (R) with finite first derivatives [cf. Eq. 

(|2.18|) ]. it follows that r is also finite and bounded away 
from zero. The first derivative of all the coordinates 
with respect to the affine parameter is always finite along 
IRNG, and therefore the spacetime is g-complete along 
IRNG. □ 

Proposition 2. Consider a spacetime with metric \2.1)) 
which obeys the Einstein equations for the stress-energy 
tensor \2.$fy . Assume that some portion of the spacetime 
undergoes complete gravitational collapse, i.e., 3t s S 
M+\{0} : R(t s ,r) = 0. If the shear-free condition fQ3|) 
is satisfied, such collapsing solutions are g-incomplete. 
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Proof. This follows from the shear-free condition 
(12.131) : for r > 0, R \ =>■ $ \ -co, which implies — by 
Eq. (|3.10jl — that f /" +oo. Consider now < e < 1, 
whence i?(i s ,e) = ee*^ = ^> $(i 8 ,e) \ -oo. By 
continuity, lim e ^o ^(^s, s) = — oo, and thus &(t s ,r) \ 
— oo, Vr € Mj. This means that IRNG are future- 
inextendible, and thus the spacetime is g-incomplete. 
We remark that, by standard extendability results for 
causal curves, one can show that future-directed timelike 
geodesies are also incomplete at R = 0. □ 

Proposition 3. A necessary and sufficient condition for 
the spacetime to contain a singularity is that the Brown- 
York quasilocal mass diverges negatively, E \ — oo. 

Proof. For r > 0, by Eq. (|2.29|) . E \ -oo implies 
$ \ — oo, and hence R \ by Eq. 12.1311 . The converse 
is also true: if R = with r > 0, then $ \ — oo and 
thus E \ -co. For r = 0, Eq. (|2.29(l and the regularity 
condition (|2.17|) imply 



lim i? = — e 



-v. 



(3.11) 



If E \ — oo then $ \ — oo, whence i? \ 0. To prove 
the converse, we first note that regularity at the cen- 
ter [Eq. (|2~T7jl ] together with Eq. (|2~T3|) imply that 
lim r ^o $ r = 0. As shown previously (cf. proof of 
Prop. 2), lim r ^o — — oo, which means that the limit 
(|3.11|) is a 0/0-undeterminacy. Successive application 
of l'Hopital's rule, together with lim r _»o = yields 
E \ — co. This completes the proof. □ 



IV. LOCAL VISIBILITY OF THE 
SINGULARITY 



Given that a spacetime singularity forms at R(t s , r) = 
0, is it (at least) locally visible? Local visibility requires 
one to show the existence of non-spacelikc future-directed 
outgoing geodesies — say, ORNG — with their past end- 
point at the singularity. Geometrically, this is equivalent 
to requiring that the area radius increases along such 
geodesies, i.e., (dRfdr)oRNG > 0. Along ORNG we have 



ar \ar 



R.r + e^Rj. (4.1) 



ORNG 



Using the standard procedure |23| , we introduce auxiliary 
variables u, X: 



u = r a , a > 0; A 



R 



(4.2) 



In the limit of approach to the singularity we have 



An 



R ,. dR 
hm — = hm — — 

R^0,u^0 u fl^o,u-+o du 



= lim 



1 dR 



R^o,r^o ar" 



lim 



1 



fl^o,r^o ar' 



dr 

(R, r + e*- 9 R t ) 



= lim 



r— >0 ar 



a-2 



4> 



H(t) 



(4.3) 



Ao is the value of the tangent to the geodesic on the 
u — X plane, at the singularity. If Ao is real and positive- 
definite, the singularity is at least locally naked, and cov- 
ered otherwise. 

Proposition 4. The R(t si 0) — singularity is also a 
shell- crossing singularity. 

Proof. Shell-crossing singularities occur when R r = 0, 
and signal the overlapping of infinitesimally close spheri- 
cal shells on a given spacelike slice E t . Such singularities 
also occur in Newtonian collapse, and, as such, are not 
intrinsic general-relativistic phenomena; in addition, for 
R > 0, they are gravitationally weak, and may be trans- 
formed away by a C 1 coordinate change (24|. To distin- 
guish shell-crossings from shell-focusing singularities, the 
condition R > is usually implied in the definition of 
the former; however, one can have R = i? r = 0, which is 
the case addressed here. Regularity at the center [cf. Eq. 
(|2.17l) ] together with the shear-free condition (|2.13|) [with 
the gauge choice f(r) = 1/r] imply lim r ^ R,r = e*^' 1 " ). 
As shown previously, lim r ^o ^K^s, r ) = —oo, and thus 
lim r _>o R,r(ts, t) = 0. □ 

Proposition 5. The R(t s ,0) = singularity in shear- 
free collapse is covered. 

Proof. The function H(t) is monotonically decreas- 
ing for spacetimes undergoing complete gravitational 
collapse [for spacetimes with a Schwarzschild exterior 
matched to the interior solution at r = r c , H(t) diverges 
negatively at t = t s (r c ), when the "last" matter shell col- 
lapses, signalling infinite convergence of the fluid flow], 
and regularity at the center requires lim r ^o = 0. 
From Eq. (|4.3|l we have then 



An 



lim 



o*(ts,r) 



+o ar a ~ 



(4.4) 



The requirement $ r — >0asr\0 implies that one 
cannot have e* oc r Q_1 . This, in turn, implies that, for 
a G (0, 1] we have Ao = 0, and the singularity is thus 
covered. For a > 1 there are two distinct possibilities 
(with the necessary assumption that e* <~k r Q_1 ): (i) 
e* > r Q_1 , or (ii) e* < r Q_1 . Possibility (i) implies that 
Ao f +oo, and (ii) yields Ao = 0; that is, for a > 1, Ao 
cannot be positive-definite, and the singularity is covered. 
□ 

Proposition 6. If the weak energy condition holds, H > 
0, and p < 1/R 2 , all noncentral singularities R(t s ,r > 
0) = are covered. 
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Proof. The sufficient condition for local visibility out- 
lined at the beginning of this section only works in the 
r \ limit, so one must resort to other methods to study 
the local visibility of nonccntral (r > 0) singularities. 
Noncentral shells will be covered if they become trapped 
before they become singular. The apparent horizon in 
spherical symmetry is given by Rit^r) = 2M(t a h, r), 
which implicitly defines the curve i a h(r). On the t — R 
plane, the condition for local coveredness is then: 



dt 



ah 



dt 



dR 



UR 
dr V dr 



ah 



all 



R, 



> o. 



(4.5) 



where t 



all 



dt a h/dr. Now, since R yt < 0, the above 



condition reads 



R,r + R,tt' ah > 0. 

Taking the total r-derivative of i?(i a h, r) 
along t a h(r) one obtains 



2M r - i? , 
R t - 2M t 



(4.6) 
(4.7) 

2M(t ah ,r) 
(4.8) 



where all quantities are evaluated along i a h(0- The con- 
dition II > together with Eq. I|2.5[l implies M t > 0, and 
p < l/R 2 together with Eq. (|2~4*j) implies 2M, r < i? r . 
This means that both the numerator and denominator 
in Eq. i|4.8|) are negative, whereby condition l|4.6[) is sat- 
isfied. Now, since R yt < and M, t > 0, condition (|4.7j) 
is equivalent to M >r R tt — M tt R, r < 0. Using Eqs. 
(|2.5|) . this reads 



R 2 R tt R,r(P + 11 ) < 0: 



(4.9) 



which holds provided the weak energy condition holds: 
T^u^u" >0^p + n>0. This completes the proof. □ 

Proposition 7. All radially null-covered R = singu- 
larities are covered along all causal directions. 

Proof. Radially null-covered R = singularities satisfy 
£ a h(f) < t*(r), where t*(r) is the solution for ORNG on 
the t — r plane, and is given by the first-order ODE: 



dt * Til \ 



(4.10) 



with the initial datum i*(0) = t s (0). Future-directed non- 
radial geodesies are given by 



dt 
dr 



e*-*4/l 



^(^+ e )>-^), (4.11) 



R 2 



with i(0) = f s (0). Since t(0) = t*(0), and (di/dr) > 
(dtsf/dr) > 0, by a straightforward application of the 
Gronwall lemma it follows that 



t(r) > U(r), Vr G Kj, 



(4.12) 



where the inequality saturates only for r = 0. Since 
^ah( r ) < t*{ r ), we also have t a h(?*) < t( r ), which means 
that the singularity is covered along all causal (radial and 
nonradial) directions. □ 



V. CONCLUDING REMARKS 

We have shown that, provided certain plausible regu- 
larity and physical reasonability conditions hold, spheri- 
cal shear-free collapse results in a locally covered singu- 
larity, which must therefore be spacelike. This may thus 
be viewed as a strong cosmic censorship result for the 
class of spacetimes considered herein. In terms of mat- 
ter content, such class includes virtually all known clas- 
sical descriptions of matter — anisotropic fluids and spe- 
cializations thereof (e.g., perfect fluids, dust), Maxwell 
fields, and scalar fields — and the constraints imposed on 
it (needed only for the particular case of noncentral shells; 
cf. Prop. 6) are arguably soft. 

What is the gravitational nature of the spacelike sin- 
gularity in shear-free collapse? In typical collapse situ- 
ations (with nonvanishing shear) leading to a curvature 
singularity, the Weyl curvature diverges faster than the 
Ricci, near the singularity (2{|. Now, shear- free collapse 
altogether forbids the influence of the Weyl tensor on 
the evolution of nonspacelike congruences, since only the 
propagation equation for shear contains a Weyl tensor 
term — the Raychaudhuri equation, which gives the evo- 
lution of geodesic expansion, is independent of the Weyl 
tensor. Since the Weyl tensor gives the purely gravi- 
tational (i.e., nonlocal) contribution, the singularity in 
shear-free collapse is essentially nongravitational in na- 
ture, being instead determined by the local matter dis- 
tribution. Indeed, order of magnitude estimates yield 
Cfi^s ~ R^v ~ r 2 for shear-free collapse, whence Weyl 
fails to dominate over Ricci [2(|. 

Finally, a natural question is: How do these results 
carry over to nonspherical geometries? The analogous 
singularity properties of spherical Tolman-Bondi collapse 
and Szekeres |22| (quasispherical) dust collapse [28j sug- 
gest that quasispherical shear-free collapse should remain 
strongly censored. For stronger departures from spheri- 
cal symmetry, given the local matter-dominated nature 
of the singularity, one may speculate that the result 
holds (see conjecture below), although this is far from 
evident. We note, however, that Szekeres spacetimes — 
notwithstanding the fact that they do not admit any 
Killing vector field — cannot contain gravitational radi- 
ation; for stronger departures from spherical symmetry, 
there are in general an infinite number of gravitational 
degrees of freedom, and the singularity might conceivably 
be Weyl-dominated. 

Based on all of the above, we conclude by putting for- 
ward the following conjecture: Let [M, g^ v ) be a space- 
time admitting an ADM decomposition M w M x S, 
where S is a submanifold with induced Riemannian 
three-metric 7^ and conjugate gravitational momentum 
iTij. Let be the stress-energy tensor in (M, g^ v ), and 
let u M be the proper four-velocity of the matter flow, 
with four-acceleration a p := u"V„!i''. If the following 
conditions hold: (i) (7ij,7ry) are regular initial data on 
S, (ii) T^i, obeys the weak energy condition, and (iii) 
C\lv ■= V^u^) - i(V Q u Q )(.g AI „ - UpUv) + a^Uyi) = 0, 
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then the Cauchy horizon in (M, g^) is the empty set. SFRH-BPD-5615-2001 and by NSF Grant PHY-0098084. 
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